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ABSTRACT 


ill  . 


The  P-codas  of  four  Aleutian  teleseismic  events,  recorded 
by  an  array  of  three  portable  seismic  stations  near  Leduc,  Alberta, 
have  been  analyzed  by  the  spectral  ratio  method  for  purposes  of 
determining  the  reliability  of  this  method  and  of  obtaining  infor¬ 
mation  about  crustal  structure. 

The  general  characteristics  of  the  spectral  ratios  were 
found  to  be  reasonably  consistent  at  each  of  the  stations  in  0.5 
to  3.0  Hz  range.  However,  a  considerable  difference  was  noted  in 
the  spectral  ratios  from  station  to  station,  a  rather  surprising 
result  since  the  stations  were  only  about  ten  kilometers  apart  in 
a  presumably  uncomplicated  geological  region  consisting  of  practi¬ 
cally  horizontal  sedimentary  layering. 

Many  theoretical  spectral  ratios  were  generated  using 
crustal  models  based  on  known  geology  along  with  varied  assumptions 
about  unknown • crustal  parameters.  While  no  exact  match  has  been 
found  at  this  time  between  theoretical  and  observed  spectral  ratios, 
it  is  possible  that  the  upper  crust  in  the  Leduc  area  is  typified 
by  1)  a  very  low  S-wave  velocity  in  the  surficial  layers,  2)  a 
Poisson's  ratio  of  one-third  in  the  softer  sedimentary  layers,  and 
3)  a  low  velocity  layer  of  thickness  approximately  3.5  kilometers 
in  the  basement  immediately  beneath  the  sediments. 

The  effect  on  the  P-coda  spectral  ratio  of  anisotropy  in 
the  sedimentary  layers  was  studied  using  a  modified  version  of  the 
Haskell— 'Thompson  matrix  method  adapted  to  transversely  isotropic 


IV 


media  for  the  case  of  body  waves  by  E.R.  Kanasewich.  Anisotropy  was 
found  to  have  a  considerable  effect  on  the  spectral  ratios  at  high 
angles  of  incidence,  but  almost  no  effect  for  low  angles  of  incidence. 
For  the  Aleutian  events  in  this  study,  sedimentary  anisotropy  was  not 
found  to  have  an  important  effect  on  the  spectral  ratios. 
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CHAPTER  1 


THE  SPECTRAL  RATIO  METHOD  OF  DETERMINING  CRUSTAL  STRUCTURE 


1.1  Introduction 

In  the  past  decade,  there  have  been  several  attempts  by  var¬ 
ious  researchers  to  use  the  spectral  characteristics  of  the  P-wave 
together  with  the  Haskell— Thompson  matrix  formulation  for  the  solution 
of  the  elastic  wave  equation  to  determine  the  nature  of  the  crust. 

Tnis  is  an  indirect  inverse  problem  in  that  It  involves  the  generation 
of  many  crustal  models  in  hopes  of  finding  a  best  fit  to  observed  data 
and  like  most  inverse  problems  in  geophysics,  the  solution  is  not  neces¬ 
sarily  unique.  Furthermore ,  the  observed  data  itself  has  been  found  to 
be  quite  susceptible  to  noise  from  sources  other  than  those  included  in 
the  Haskell-Thompson  theory.  Because  of  these  problems,  the  degree  to 
which  the  Haskell-Thompson  formulation  is  capable  of  resolving  the 
crustal  structure  still  remains  undetermined.  The  purpose  of  this 
thesis  is  to  find  to  what  extent  the  application  of  the  Haskell- 
Thompson  matrix  to  crustal  structure  is  possible  and,  then,  to  actually 
apply  the  method  to  determine  the  crustal  structure  as  accurately  as 
possible  in  central  Alberta. 

When  the  crust  of  the  earth  is  penetrated  by  seismic  waves , 
its  frequency  response  to  the  signal  depends  on  the  elastic  parameters 
and  thicknesses  of  the  rock  layers  within  the  crust.  If  the  crust  is 
assumed  to  be  composed  of  strictly  horizontal  layers,  a  condition  gener¬ 
ally  verified  by  refraction,  surface  wave  and  teleseismic  travel  time 
studies,  the  frequency  response  or  transfer  function  of  the  earth  to 
plane  P  or  S  waves  can  be  computed  by  solving  the  wave  equation  using 
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the  continuity  of  stresses  and  particle  displacements  at  each  interface 


and  the  vanishing  of  stresses  at  the  free  surface  as  boundary  conditions 
The  I  la sk el 1 -Thompson  matrix  formulation  performs  this  calculation  for 
any  number  of  horizontal  layers  in  terms  of  a  A  x  A  matrix  whose  members 
are  a  function  of  the  parameter  of  each  layer  and  the  angle  of  incidence 
of  the  source  of  excitation. 

The  crust  acts  like  a  linear  filter  on  incoming  seismic  waves 
and  its  effect  on  resultant  seismograms  is  quite  significant  and  has 
been  observed  at  seismological  stations  for  many  years.  Gutenberg 
(1934)  found  that  certain  seismic  stations  showed  a  tendency  to  detect 
seismic  waves  of  a  given  period.  Neumann  (1954)  noticed  that  stations 
located  on  deep  sedimentary  rocks  recorded  less  amplitude  for  a  given 
earthquake  than  those  stations  located  directly  on  the  granite  basement. 
The  effect  of  the  upper  crust  on  the  amplitudes  of  body  waves  led 
Gutenberg  and  Richter  (193G)  to  devise  station  correction  factors  for 
earthquake  magnitude  determinations. 

It  has  been  only  recently,  however,  that  attempts  have  been 
made  to  invert  the  observed  crustal  effects  on  seismic  body  waves  and, 
hence,  arrive  at  the  crustal  structure.  Such  a  method  would  be  quite 
valuable  to  crustal  studies  in  that  the  only  other  ways  to  detail  the 
structure  of  the  crust,  namely,  refraction  and  reflection  seismology, 
are  plagued  with  the  necessity  of  artificial  sources,  extensive  field 
operations,  and  accurate  time  correlations  between  stations.  An  inver¬ 
sion  method,  however,  based  on  the  application  of  the  Haskell-Thompson 
matrix  formulation  to  earthquake  body  waves  would  be  free  from  the  above 
mentioned  encumber ances  and  would ,  thus ,  be  a  very  inexpensive  method 


of  obtaining  crustal  information.  Phinney  (1964)  was  the  first  to  use 
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the  method  when  he  compared  the  ratio  of  long  period  spectra  of  vertical 
and  horizontal  P-codas  of  large  earthquakes  with  the  ratio  of  vertical 
and  horizontal  transfer  functions  as  calculated  by  the  Haskell-Thompson 
matrix  formulation.  His  results  for  such  long  periods  (up  to  0.2  Hz) 


wen 


quite  reasonable.  Fernandez  (1966)  also  obtained  good  fits  between 


observed  and  calculated  spectral  ratios  in  the  same  frequency  range  for 
simple  2-layer  models  (one  layer  and  half  space).  He  noted  that  any 
estimation  of  different  layers  within  the  crust  was  neither  "precise 
nor  unique"  and  suggested  that  spectra  of  shorter  period  records  would 
give  much  more  detailed  information  about  the  structure  within  the 
crust.  Ellis  and  Basham  (1968)  attempted  the  inversion  for  shorter 
period  (0.1  to  3.0  Hz)  data  in  the  central  Alberta  plains  where  the 
crustal  structure  was  presumably  well  known  but  obtained  disappointing 
agreement  between  the  theoretical  and  observed  spectra.  Alpaslan  (1968) 
also  used  short  period  data  in  such  an  inversion  process  but  without 
much  success.  Hasegawa  (1971)  obtained  fairly  promising  results  for 
the  method  at  both  long  and  short  period  results  over  a  relatively 
simple  crust  at  Yellowknife,  N.W.T. 

The  University  of  Alberta,  at  the  present  time,  has  in  oper¬ 
ation  nine  portable  broad  band  seismic  systems  capable  of  recording 
both  long  and  short  period  earthquake  waves  simultaneously.  It  is 
anticipated  that  this  system  will  be  useful  in  determining  the  actual 
limitations  and  possibilities  in  inverting  P-wave  spectral  data. 


1.2  The  Inversion  of  the  Spectral  Ratio 

A  signal  transmitted  through  the  crustal  layers  of  the  earth 
and  recorded  at  the  surface  is  modified  in  a  manner  depending  on  the 


density  and  elastic  properties  of  the  media. 
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Such  a  recorded  signal  is 
expressible  in  the  time  domain  as  the  convolution  of  the  source  function 
with  the  crustal  transfer  function.  That  is, 

Recorded  Seismograms  (t)  =  Source  (t)  *  Transfer  Function  (t)  1.1 

where  *  indicates  a  convolution  process  and  t  represents  time.  It  is 
far  more  convenient  to  consider  the  above  relation  in  the  frequency 
domain  for  several  reasons:  1)  it  is  easy  to  calculate  the  crustal 
transfer  as  a  function  of  frequency  using  normal  mode  theory  via  the 
Haskell— Thompson  matrix  formulation,  2)  statistical  methods  of  reducing 
recorded  ground  motion  to  its  spectral  components  are  readily  available, 
and  3)  convolution  in  the  time  domain  becomes  simple  multiplication  in 
the  frequency  domain .  In  the  frequency  domain ,  therefore : 

Spectrum  of  recorded  data  =  Source  *  Crustal  Transfer  Function  1.2 

It  should  be  noted  that  the  term  spectrum  includes  both  phase  and 
amplitude  Information  as  a  function  of  frequency,  hence,  the  crustal 
transfer  function  itself  Is  a  complex  quantity  for  each  frequency. 

If  one  is  considering  only  plane  P~waves  as  the  source,  all  of 
which  are  incident  on  the  base  of  the  crust  at  the  same  angle,  then  the 
resultant  motion  on  the  surface  will  be  only  in  the  plane  determined  by 
the  radial  and  vertical  directions.  Hence  the  recorded  data  can  be 
specified  In  terms  of  two  components  only:  the  vertical  and  horizontal 
motion. 


Since  both  the  vertical  and  horizontal  ground  motion  are  due 


to  exactly  the  same  source  one  can  write: 

Z  (f)  =  Source  (f)  *  V  (f) 

H  (f)  =  Source  (f)  -  R  (f) 

where :  Z  and  H  are  the  vertical  and  horizontal  spectrums  of  the  r 

corded  data;  V  and  R  are  the  vertical  and  horizontal  crustal  trass 
functions;  and  f  represents  frequency.  Simple  division  of  these  f 
equations  leads  to  the  basis  for  the  inversion  of  P-wave  spectral 

Z  (f)  V  (f) 

H  (f)  R  (f) 

that  is ,  the  ratio  of  the  vertical  to  horizontal  spectrums  of  the 
'corded  seismograms  is  equal  to  the  ratios  of  the  vertical  to  horiz 
transfer  functions.  The  source  function  is,  therefore,  eliminated 
the  division  process,  and  the  observed  spectral  ratio  can  be  match 
the  theoretically  calculated  ratio  of  transfer  functions  for  any  a 
crustal  model. 

In  practice,  however,  several  difficulties  arise.  In  re 
to  the  eliminations  of  the  source,  one  must  recall  the  assumptions 

it: 

1)  It  must  be  a  plane  wave.  From  a  geometric  point  of  v 
this  assumption  should  be  quite  valid  for  most  teleseismic  earthqu 
but  any  large  lateral  heterogeneities  in  the  mantle  below  the  crus 
could  alter  the  plane  characteristics  of  the  wave,  i.e.  an  active 


relict  subduction  zone. 
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2)  The  source  of  plane  waves  impinging  on  the  layered  media 
must  consist  of  pure  compressional  (or  shear)  motion.  Because  of  the 
difference  in  P  and  S  velocities,  the  distance  between  the  source  origin 
and  the  point  of  observation  acts  as  a  delay  line  to  separate  P  and  S 
motion.  Standard  travel  time  curves  are,  therefore,  useful  for  deter¬ 
mining  how  much  of  the  recorded  data  at  particular  epicentral  distance 
should  be  P  motion.  The  real  danger  of  S  waves  entering  the  P-coda 
comes  not  from  shear  waves  generated  near  the  source  but  rather  from 
shear  waves  produced  by  P-S  conversions  in  the  mantle  beneath  the  crust. 
While  there  is  no  seismic  evidence  of  significant  first  order  velocity 
jumps  in  the  mantle,  there  are  large  velocity  increases  especially  in 
the  100  to  700  KM  depth  range  which  might  look  like  discontinuities  to 
seismic  waves  of  certain  wavelengths . 

3)  All  the  components  of  the  source  wave  must  arrive  at  the 
base  of  the  crust  at  the  same  angle  of  incidence.  This  assumption 
eliminates  the  use  of  relatively  close  events  as  sources  since  it  is 
quite  likely  that  P-waves  recorded  at  the  station  have  arrived  by  many 
different  paths.  Aipaslan  (1968)  attempted  to  use  Project  Early  Rise 
explosion  data  to  invert  the  P-wave  spectral  ratios,  but  did  not  have 
much  success  partly  because  of  the  varying  angles  of  incidence  of  the 
incoming  P-waves  as  they  were  refracted  or  reflected  along  different 
mantle  boundaries.  It  is,  therefore,  important  that  the  events  chosen 
for  analysis  be  at  a  safe  distance  from  the  recording  station  to  avoid 
these  multiple  arrivals, 

4)  Another  possible  cause  for  the  angle  of  incidence  not 
being  constant  at  the  base  of  the  crust  for  the  entire  P-coda  is  the 
possible  presence  of  first  order  discontinuities  in  the  mantle.  Any 
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multiple  reflection  in  the  mantle  will  affect  some  frequencies  more 
than  others,  and  the  angle  of  emergence  of  the  source  wave  from  the 
mantle  which  impinges  on  the  crust  would  be  a  function  of  frequency, 
rather  than  a  constant  as  is  assumed  in  the  theory.  this  effect  can 
be  reduced  by  smoothing  of  the  spectral  estimates  because  the  large 
time  lags  involved  when  considering  mantle  structure  appear  as  a  high 
frequency  modulation  on  the  spectra. 

1.3  Observed  Spectral  Ratios 

The  inversion  of  P-wave  spectral  data  involves  the  calculation 
of  the  ratio  of  the  vertical  and  horizontal  components  of  ground  motion 
in  the  frequency  domain.  This  ratio  involves  both  phase  and  amplitude 
information  and  is,  therefore,  a  complex  quantity;  so  for  comparison 
purposes,  the  modulus  or  absolute  value  of  the  ratio  is  generally 
considered  alone. 

A  simple  Fourier  transform  of  the  seismogram  would  yield  its 
complex  spectrum,  however,  for  reasons  considered  below,  it  is  more 
convenient  to  consider  the  power  spectrum  of  the  seismic  data.  As  noted 
above,  the  mantle  of  the  earth  is  possibly  a  cause  of  multiple  reflec¬ 
tions  and  conversions  of  teleseisrnic  wave  prior  to  their  arrival  at  the 
base  of  the  crust.  Such  a  mantle  structure  has  large  time  lags  associ¬ 
ated  with  it,  and  close  spectral  peaks,  while  the  effect  of  crustal 
layers  on  the  spectrum  should  be  small  time  lags  and  widely  spaced 
spectral  peaks.  Figure  A-l  shows  theoretical  spectral  ratios  for  simple 
2-layer  models  of  various  thicknesses,  but  with  identical  layer  veloc¬ 
ities.  The  thicker  the  model,  the  closer  the  spectral  peaks  approach 
each  other.  The  use  of  the  power  spectrum  in  the  analysis  rather  than 
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the  simple  Fourier  transform  allows  one  to  smooth  the  resulting  spectrum 
and  thus  eliminate  any  mantle  structure  from  the  problem.  Ellis  and 
Basham  (1968)  carried  this  smoothing  procedure  even  farther  and  pro¬ 
ceeded  to  remove  not  only  the  effect  of  mantle  structure  from  the  spec¬ 
trum  but  also  the  effect  of  the  whole  lower  crust. 

In  practice,  the  maximum  time  lag  which  can  be  considered  in 
the  power  spectrum,  is  equal  to  the  length  of  the  P-coda  under  consid¬ 
eration  because  the  power  spectrum  is  a  Fourier  cosine  transform  of 
the  autocorrelation  function  of  the  data.  Only  lags  from  zero  to  some 
finite  value  are  available  so  the  autocorrelation  is  only  an  estimate. 
One  must,  therefore,  effectively  multiply  the  true  autocorrelation  with 
infinite  lag  by  a  time  windovr  of  the  same  length  as  the  available  data. 

1.4  Theoretical  Spectral  Ratios 

Since  a  time  window  must  be  applied  to  the  observed  data,  a 
similar  window  should  be  applied  to  the  theoretically  determined  spec¬ 
tral  ratios.  LeBlanc  (1966),  the  first  to  study  this  problem  in  detail, 
showed  that  comparing  the  theoretical  truncated  spectral  ratio  to  the 
observed  truncated  spectral  ratio  was  a  valid  process  provided  the 
source  was  short  compared  to  the  record  length.  Alpaslan  (1968)  gave  a 
straight  forward  means  of  generating  truncated  theoretical  spectral 
ratios.  He  convolved  an  assumed  source  with  the  theoretical  transfer 
functions  to  yield  vertical  and  horizontal  component  synthetic  seismo¬ 
grams.  These  synthetics  were  then  treated  precisely  in  the  same  manner 
as  the  observed  seismograms,  that  is,  the  same  time  window  was  chosen, 
and  the  power  spectral  ratios  calculated  in  the  same  way.  This  process 
results  In  the  comparison  of  theoretical  and  experimental  results  both 


of  which  have  undergone  exactly  the  same  type  of  processing.  This 
method  was  used  in  this  study  although,  for  the  length  of  window  used 
(forty  seconds)  and  in  the  frequency  range  of  0.5  to  4.0  Hz,  the 
difference  between  the  truncated  and  non-truncated  theoretical  spectral 
ratio  was  found  to  be  actually  insignificant. 
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CHAPTER  II 
EXPERIMENTAL  DATA 


2.1  Introduction 

In  the  fall  of  1969  an  array  of  three  seismograph  stations 
was  operated  by  The  University  of  Alberta  in  the  area  immediately  south¬ 
east  of  the  town  of  Leduc,  Alberta.  Each  station  consisted  of  a  three 
component  broad  band  seismic  system  recording  directly  on  seven  track 
digital  tape.  The  system  is  described  in  detail  in  Burke,  Kanasewich, 
Malinsky,  and  Montalbetti  (1970).  Basically  it  is  composed  of  Willmore 
MK  II  seismometers  with  natural  periods  of  1.5  seconds  and  an  electronic 
amplifier  which  extend  the  usable  velocity  sensitivity  out  to  the  40 
second  period  range  under  normal  conditions.  This  system  response  is 
obtained  by  the  use  of  second  order  active  RC  network  filters  and  1 ow 
noise  amplifiers  which  separate  long  and  short  period  data,  amplify  the 
long  period  data  to  compensate  for  the  drop  in  sensitivity  due  to  the 
seismometer,  and,  then,  recombine  the  two  signals  in  such  a  way  as 
either  to  enhance  or  attenuate  microseismic  noise.  An  analog  to  digital 
converter  was  used  to  record  the  data  on  a  seven  track  digital  field 
tape. 

The  amplifiers  were  calibrated  and  matched  in  the  laboratory 
but  no  calibration  was  undertaken  of  the  whole  system  including  the 
seismometers.  Hence  there  is  some  uncertainty  as  to  the  relative 
amplitudes  of  the  components  of  the  various  events  at  the  three  field 
stations.  Careful  calibration  must  certainly  be  done  in  the  future  for 
spectral  ratio  work  which  involves  amplitude  ratios  of  different  systems 


at  different  frequencies. 
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The  microse ramie  attenuation  feature  of  the  system  proved  to 
be  unfortunate  from  the  point  of  view  of  studying  spectral  ratios  and 
no  useful  long  period  spectral  ratio  data  greater  than  two  second 
period  was  obtained  in  the  fall  of  1959.  The  higher  frequency  range 
(0.5  to  3.0  Hz)  produced  good  data,  however,  particularly  for  studies 
of  the  upper  crust. 

Four  events,  all  from  the  Aleutian  Islands,  were  chosen  for 
spectral  ratio  studies.  Pertinent  information  about  these  events  is 
shown  in  Table  2-1.  Inspection  of  the  seismograms  show  them  all  to  have 
a  high  signal  to  noise  ratio.  Since  all  the  events  come  from  the  same 
area,  no  discrepancies  from  either  difference  in  angles  of  incidence 
or  azimuthal  inhomogeneities  around  the  station  are  expected.  Inspection 
of  travel  time  charts  shews  that  for  the  epicentral  distance  of  about 
forty  degrees,  no  S-waves  generated  at  the  source  can  be  expected  for 
almost  six  minutes  after  the  onset  of  the  P-waves;  the  PP-wave  is 
delayed  about  one  hundred  seconds;  and  the  PcP  over  two  minutes.  Hence 


these  phases  are  not  expected  to  interfere  with  the  initial  P-coda. 


TABLE  2-1 


EVENT 

DATE 

Amchitka 

2  Oct. 

Nuclear- 

1969 

Test  No.  1 

A1 eutians 

18  Oct 

OvJ 

• 

0 

y~zr 

1969 

Andreanof f 

21  Oct 

No.  3 

1969 

Aleutians 

31  Oct 

No.  4 

1969 

MAGNITUDE  EPICENTRAL  DEPTH 
DISTANCE 

6.5  39.0°  - 


5.6  42.3°  24  KM 

5.9  39.1°  48  KM 


39.0 


o 


AZIMUTH  STATIONS 
AT  OBS  RECORDING 


294°  OBS,  VNG 

295°  OBS,  OSC, 

VNG 

299°  OBS,  OSC, 

VNG 

294°  OSC,  VNG 


6.0 


49  KM 
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2.2  'The  Effect  of  the  Phase  pP 

The  earthquakes  chosen  occur  at  depth,  and  pP  phases  are 
expected  only  seconds  after  the  P  impulses.  Phinney  (1964)  has  studied 
the  effect  of  pP  phases  on  the  spectral  ratios  and  concluded  that, 
provided  the  amplitude  of  these  phases  are  small  in  comparison  with  the 
initial  pulse,  they  will  not  affect  the  validity  of  the  analysis  since 
the  signal  will  retain  its  random  nature.  If  the  pP  phase  is  signif¬ 
icantly  large,  however,  a  perturbation  ripple  with  a  frequency  spacing 
of  one  over  the  pP-P  time  lag  will  show  up  on  the  spectrum.  To  test 
for  the  existance  of  this  perturbation  on  the  spectral  ratios ,  an 
experiment  was  performed  on  data  obtained  at  the  OSC  station.  A  cepstrum 
analysis  of  the  observed  spectral  ratio,  as  used  by  Bakun  (1971),  was 
performed  to  separate  those  recurring  features  having  a  close  frequency 
spacing  from  those  having  a  longer  frequency  spacing.  It  is  interesting 
to  note  that  such  an  analysis  has  units  of  cycles  per  hertz  or,  simply, 
seconds.  Hence,  one  can  obtain  from  the  cepstrum  a  graph  of  the  ampli¬ 
tude  of  the  spectral  ratio  versus  the  time  lags  which  have  presumably 
produced  the  characteristics  of  the  transfer  function.  Short  time  lags 
would,  for  instance,  arise  from  the  sedimentary  layers  of  the  crust, 
longer  time  lags  from  the  structure  of  the  deep  crust,  and  the  longest 
time  lags  from  the  mantle  structure.  If  the  pP  phase  significantly 
perturbs  the  spectral  ratio,  then,  much  information  from  the  lower  part 
of  the  crust  might  be  lost  due  to  interference  from  the  pP  effect. 

Figure  A— 2  shows  the  cepstrum  curve  for  the  OSC  events.  Two  of  the 
events,  OSC  3  and  OSC  4  occurred  at.  depths  of  48  KM,  thus,  giving  rise 
to  a  pP  phase  arrival  at  about  twelve  seconds  after  the  initial  P. 

This  phase  can  be  identified  on  the  records  quite  well.  The  curve  does 
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show  peaks  for  the  OSC  3  and  OSC  4  events  at  the  predicted  time  lag  of 
about  twelve  seconds.  The  OSC  2  event,  however,  which  occurred  at  a 
depth  of  only  24  KK  with  a  pP-P  lag  of  six  seconds. has  a  quite  different 
curve  for  time  lags  greater  than  five  seconds  which  suggests  that  the 
similarity  between  the  curves  of  OSC  3  and  OSC  4  is  not  due  to  the 
similar  lower  crust  beneath  them  since  the  crustal  effect  would  have 
influenced  OSC  2  as  well.  The  pP-P  effect  on  the  spectral  ratios  is, 
thus,  seen  to  be  significant,  and  may,  indeed,  be  a  more  important 
source  of  the  close-spaced  ripples  on  the  ratios  than  the  lower  crustal 
structure.  It  would  be  possible  to  avoid  this  effect  by  choosing  for 
spectral  ratio  studies  only  those  events  which  are  very  shallow  such  as 
nuclear  tests  or  those  events  which  are  very  deep  such  that  the  pP-P 
time  lag  would  be  so  large  as  to  not  affect  the  P-coda  window  at  all. 

In  this  study,  concerned  mostly  with  only  the  effect  of  the  sediments 
on  the  spectral  ratio,  the  perturbation  ripples  due  to  the  pP-P  lag  are 
no  larger  than  0.2  Hz  which  would  not  affect  the  general  characteristics 
of  the  spectral  ratios.  In  fact,  Figure  A- 2  shows  considerable  simi¬ 
larity  between  all  three  events  for  time  lags  less  than  five  seconds 
which  would  be  sufficient  for  any  effects  of  the  upper  crustal  layers. 

2.3  Rotation  of  the  Observed  Data 

Since  the  data  v/as  recorded  using  a  tri-partite  system  aligned 
in  the  vertical,  east-west,  and  north-south  directions,  the  horizontal 
components  of  ground  motion  were  first  rotated  to  obtain  a  radial  and 
transverse  record.  To  accomplish  this,  the  expected  great  circle 
azimuth  from  source  to  receiver  v/as  computed  and  a  transformation  was 
applied  to  the  data  as  follows : 
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R  -  E  sin  0  -i-  N  cos  0 


-E  cos  0  -i-  N  sin  0 


2.1 


where  0  is  the  great  circle  azimuth  measured  clockwise  from  the  north, 


R,  T,  E,  and  N  are  the  radial,  transverse,  oast-west,  and  north-south 
components  of  ground  motion  respectively. 

Ideally,  the  P-coda  should  produce  ground  motion  in  the 
vertical  and  radial  directions  only  for  a  crust  consisting  of  horizontal 
layers.  A  certain  degree  of  scattering  in  the  crust,  however,  must  be 
expected  and,  hence,  there  is  motion  in  the  transverse  plane  on  the 
observed  records.  The  transverse  motion  is  generally  small  compared 
to  the  vertical  and  radial  motion,  however,  and  hopefully,  the  scat¬ 
tering  does  not  affect  the  spectral  ratios  drastically. 

Ellis  and  Basham  (1963),  by  means  of  particle  motion  studies, 
have  suggested  that  seismic  waves  arriving  at  Leduc  are  coming  in  at 
an  azimuth  averaging  18  degrees  further  north  than  predicted  from  the 
great  circle  path.  An  experiment  was  performed  to  test  this  assertion 
by  an  independent  method  and  also  to  obtain  more  information  about  the 
transverse  motion  on  the  records.  The  north  and  south  components  were 
rotated  in  one  degree  increments  and  the  maximum  value  of  the  auto¬ 
correlation  function  of  the  transverse  trace  was  calculated  for  each 
azimuth  angle.  Since  the  autocorrelation  function  gives  an  estimate 
of  the  amount  of  energy  in  a  trace,  that  azimuth  angle  for  which  the 
autocorrelation  is  minimum  yields  an  apparent  azimuth  angle  for  each 
event.  Table  2-2  shows  those  azimuth  angles  for  which  the  maximum 


value  of  the  autocorrelation  function  was  a  minimum  for  the  various 


events  at  the  three  field  station 
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TABLE  2-2 


EVENT  STATION  ANGLE  OF  MINIMUM 

TRANSVERSE  ENERGY 


Amchitka 

VNG 

318° 

Aleutians 

VNG 

337° 

18  Oct. 

OBS 

255° 

OSC 

295° 

Andreanof f 

VNG 

268° 

OBS 

264° 

OSC 

306° 

Aleutians 

VNG 

293° 

31  Oct. 

OSC 

293° 

The  expected  great  circle  azimuths  for  these  events  are  in 
the  295°  to  300°  range.  It  is  likely  that  since  only  a  hand  compass 
was  used  to  orient  the  stations  at  OSC  and  VNG,  there  might  be  several 
degrees  error  inherent  in  these  results,  but  the  amount  of  scatter 
indicates  that  this  method  is  probably  not  a  very  reliable  means  of 
obtaining  the  actual  azimuth  angle.  However,  it  is  interesting  to 
note  the  very  low  azimuth  angles  at  OBS.  This  phenomenon  is  perhaps 
indicative  of  a  large  amount  of  transverse  energy  being  generated  in 
the  immediate  locality  of  this  station. 

2.4  Spectral  Analysis  of  the  Field  Data 

Analysis  of  the  rotated  records  indicated  that  useful  P-coda 
information  was  arriving  at  least  forty  seconds  after  the  first  P  pulse. 
At  times  later  than  forty  seconds,  the  transverse  trace  appeared  to 
have  as  much  energy  as  the  vertical  and  radial  components  and  it  was 
not  felt  that  this  later  information  would  be  very  useful.  Hence,  the 
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first  forty  seconds  of  each  event  were  used  to  obtain  spectral  ratios. 

To  avoid  a  sharp  edge  on  the  end  of  the  data,  a  cosine  taper  was  used 
on  the  last  four  seconds  of  the  forty  second  window.  Actually,  exper¬ 
iments  showed  that  for  a  window  this  long  the  taper  had  very  little 
effect  on  the  final  spectral  ratios.  Since  a  forty  second  P-coda 
necessarily  includes  the  pP  phase  for  medium  depth  Aleutian  events,  the 
spectral  ratios  were  also  calculated  for  a  ten  second  P-coda.  The  ten 
second  window,  however,  provided  poor  resolution  of  the  spectral  ratio. 
Figure  A-3  shows  the  comparison  for  the  October  31  event  of  spectral 
ratios  with  a  ten  second  window  and  a  forty  second  window.  Since  some 
of  the  larger  peaks  seen  with  the  forty  second  v/indow  are  quite  narrow, 
it  is  reasonable  to  expect  they  might  be  diminished  in  the  low  resolution 
involved  when  one  considers  only  ten  seconds  of  data  digitized  at  0.08 
second  intervals. 

The  forty  second  P-coda  contained  501  points  to  which  an  addi¬ 
tional  523  zeros  were  added  to  bring  the  total  number  of  points  in  the 
time  sequence  to  1024  for  each  trace.  The  autocorrelation  of  the  ver¬ 
tical  and  horizontal  components  were  computed.  A  Parzen  window  was 
applied  to  the  autocorrelations  for  truncation  purposes  and,  finally, 
the  inverse  Fourier  transform  of  these  modified  autocovariances  was 
performed  to  yield  vertical  and  horizontal  power  spectrums  for  512 
points  ranging  from  .01  Hz  to  6.25  Hz.  Fourier  transforms  and  their 
inverses  were  computed  using  a  high-speed  algorithm  developed  by 
Gentleman  and  Sande  (1966).  The  square  root  of  the  ratio  of  the  ver¬ 
tical  to  horizontal  spectrums  at  each  frequency  increment  was  then 
computed  to  yield  the  observed  spectral  ratios.  Since  it  is  possible 
that  a  particular  event  may  have  no  energy  at  certain  frequencies,  a 


rejection  criteria  should  be  used  tc  avoid  these  frequencies.  For  the 
few  events  used  in  this  study,  both  the  vertical  and  horizontal  spec- 
trums  were  studied  before  the  ratios  were  computed  and  no  frequency 
voids  were  apparent. 

2.5  The  Reliability  of  the  Observed  Data 

An  interesting  fact  about  spectral  ratio  studies  is  that  at 
a  particular  station,  .totally  different  events  at  approximately  equal 
epicentral  distances  should  produce  exactly  the  same  P-coda  spectral 
ratios  regardless  of  the  nature  of  the  individual  earthquakes.  The 
actual  observed  spectral  ratios  for  the  aforementioned  events  at  the 
three  different  stations  are  shown  in  Figures  A-4,  A-5,  and  A-6.  The 
extent  to  which  the  spectral  ratios  are  alike  at  each  station  is  a 
measure  of  the  reliability  of  the  spectral  ratio  method.  For  the 
fall  1969  data,  it  is  clear  that  at  least  the  general  character  of  the 
spectral  ratios  are  significant  at  each  of  the  three  stations.  The 
smaller  "second  order"  peaks  and  troughs  are  less  useful  in  that  they 
change  for  the  different  events. 

The  reliability  of  the  Leduc  Observatory  (OBS)  data  can  be 
checked  with  data  from  other  sources.  Ellis  and  Basham  (1968)  and 
Jensen  (1970)  have  previously  recorded  events  at  the  Leduc  Observatory 
in  the  forty  degree  epicentral  range  and  have  calculated  the  spectral 
ratios.  Figure  A- 7  shows  their  spectral  ratios  compared  with  the 
Leduc  ratios  from  this  study.  The  agreement  is  quite  good  considering 
the  large  amount  of  smoothing  used  by  the  above  researchers  in  their 


analyses. 


^inc a  it  is  desirable  to  hove 


single  observed  spectral  rati 


o 


for  each  station  for  interpretat tonal  purposes,  a  non- linear  filter  v;as 
designed  to  combine  the  spectral  ratios  from  the  different  events 
according  to  the  follov.dnq  formulae: 
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2.2 


The  effect  of  this  method  of  combining  the  data  from  the  various  events 
is  to  remove  from  the  combined  data  those  peaks  and  troughs  which  are 


not  common  to  all  the  stations. 


The  stronger  the  combination,  the  more 
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drastic  is  the  removal  of  dissimilar  character.  For 
of  one,  the  combination  reduces  to  simply  an  average 
combined. 


a  filter  strength 
of  the  events 


CHAPTER  III 


THE  CRUSTAL  MODEL  FOR  CENTRAL  ALBERTA 


3.1  Introduction 


In  order  to  generate  a  theoretical  spectral  ratio  for  ar. 
given  crustal  model,  it  is  necessary  to  know  the  thickness,  P  and 
velocities,  and  density  in  each  assumed  layer  as  well  as  the  angle 
incidence  of  the  source  wave  at  the  base  of  the  bottom  layer.  The 
complexity  of  the  crust  in  central  Alberta  makes  the  determination 
of  many  of  these  parameters  difficult.  Even  the  basic  number  of  1 
required  to  properly  represent  the  crust  is  not  really  known.  Ihi 
chapter  is  an  attempt  to  resolve  some  of  these  difficulties  and  ar 
at  a  suitable  crustal  model  for  central  Alberta. 

This  region  is  geologically  characterized  by  a  thick  fla 
bedded  sedimentary  section  dipping  generally  less  than  one  hundred 
feet  per  mile  to  the  southwest.  In  the  area  of  The  University  of 
Alberta’s  Leduc  Observatory,  (Figure  A-8),  the  lithology  of  the  up 
2  KM  of  the  sediments  is  well  known  due  to  extensive  exploratory 
drilling  into  the  Devonian  reef  zones.  The  Paleozoic  section  bene 
the  reef  zone  is  much  less  known  since  few  wells  are  drilled  into 
non-oil  bearing  depths.  Three  wells  in  the  Edmonton  area  have  pen 
the  Precambrian  basement.  Interpolation  between  these  wells  gives 
total  thickness  of  the  sediments  at  the  Leduc  Observatory  to  be  at 
2.7  KM,  as  shown  by  the  basement  structure  contours  in  Figure  A— 3. 
Only  two  wells  in  the  immediate  vicinity  of  the  observatory'  have  h 
sonic  logs  run  in  them.  One  of  these  wells  is  fortuitously  locate 
almost  on  the  line  between  the  Leduc  Observatory  (0B5)  and  the  two 
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temporary  auxiliary  stations  OSC  and  VNG,  hence,  good  control  on  P-v/ave 
velocity  is  expected.  Unfortunately  the  sonic  logs  do  not  cover  the 
upper  200  meters  of  the  section  nor  any  beds  below  the  reef  zone,  so 
additional  velocity  information  was  obtained  from  reflection  seismograms 
in  the  area.  No  published  information  is  available  for  S-wave  veloc¬ 
ities,  depth  of  weathering  or  layer  densities. 

Interpolation  of  the  well  data  in  the  area  yields  the  struc¬ 
ture  through  the  three  stations  as  tabulated  in  Table  3-1.  The  P-v/ave 
velocities  down  to  the  Cooking  Lake  formation  are  an  average  of  the  two 
available  sonic  logs  with  more  weight  being  given  to  the  log  nearest 
the  stations.  The  P  velocities  In  the  lower  sediments  are  from  a 
seismic  reflection  survey  shot  in  the  area.  The  basement  velocity  is 
a  value  obtained  from  a  regional  refraction  survey  in  southern  Alberta 
(Kanasewich  and  dimming,  1966). 

Formation  density  information  is  not  readily  available.  The 
effect  of  density  on  the  spectral  ratios  is  strictly  one  of  amplitude, 
since  the  density  affects  the  reflection  coefficients  of  the  various 
boundaries  but  not  the  speed  nor  direction  with  which  the  waves  travel 
through  the  layers.  Figure  A- 9  shows  two  theoretical  ratios  for  typical 
nine  layer  Alberta  sedimentary  crustal  models,  the  only  difference 
between  the  models  being  the  density  distribution.  The  different 
density  models  produce  identical  ratios  except  for  amplitude  variations. 
Because  density  does  not  alter  the  character  of  the  theoretical  ratios , 
the  P-v/ave  velocity-density  curve  of  Nafe  and  Drake  was  used  to  assign 
density  values  to  the  crustal  layers. 

The  S-wave  velocities  in  the  sediments  of  central  Alberta  are 
not  known.  The  visual  assumption  of  rocks  behaving  as  Ideal  Poisson 
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TABLE  3-1 


TH 

ICKNESS 

( KH ) 

P— WAVE  ^ 

velocity" 

DENSITY 

(GM/CC) 

OBS 

03  C 

VNG 

(KM/ SEC) 

.01 

.01 

.01 

2.1 

2.1 

.60 

.62 

.64 

2.3 

2.1 

.19 

.21 

.21 

2.9 

2.2 

.10 

.11 

.11 

3.15 

2.2 

.04 

.04 

.04 

2.85 

2.2 

.04 

.04 

.04 

3.30 

2.3 

.08 

.08 

.08 

2.85 

2.2 

.04 

.04 

.04 

3.70 

2.3 

.02 

.02 

.02 

2.80 

2.2 

.13 

.14 

.16 

3.40 

2.3 

.07 

.08 

.08 

3.85 

2.3 

.15 

.15 

.19 

5.80 

2.6 

.24 

.18 

.24 

4.00 

2.4 

.21 

.21 

.21 

5.20 

2.6 

.21 

.21 

.21 

5.00 

2.6 

.52 

.56 

.58 

4.30 

2.5 

_ _ 

_ . _ , 

6.10 

2.7 

TABLE  3-2 

NAME 

THICKNESS 

(KM) 

P— WAVE 

VELOCITY 

(KM/SEC) 

DENSITY 

(GR./CC) 

POST  ALBERTA 

.61 

2.30 

2.1 

ALBERTA 

.51 

3.00 

2.2 

BLAIRI-iORE 

.20 

3.54 

2.3 

WABAMUN 

.15 

5.80 

2.6 

IRETON 

.24 

4.00 

2.4 

CK  LK 

.21 

5.20 

2.6 

ELK  PT 

.19 

5.00 

2.6 

CAMBRIAN 

.52 

4.30 

2.5 

BASEMENT 

- - - 

6.10 

2.7 
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solids  with  the  P-wave  velocity  equal  to  1.73  of  the  S-wave  velocity 
was  not  necessarily  considered  reliable.  Trie  problem  of  assigning 
S-wave  velocities  will  be  considered  in  detail  later.  As  a  first 
approximation,  the  S-waves  are  given  a  velocity  equal  to  one  half  of 
that  of  the  P-waves.  This  ratio  corresponds  to  a  solid  having  a 
Poisson’s  ratio  of  one- third. 

The  initial  sedimentary  model  involving  seventeen  layers  is 
based  on  the  analysis  of  the  sonic  logs.  Because  of  the  thinness  of 
some  of  these  layers,  it  was  felt  they  would  have  little  effect  on 
the  spectral  ratios  in  the  frequency  range  in  question  (0.2  -  4.0  Hz), 
so  the  seventeen  layer  crust  was  reduced  to  a  nine  layer  crust  as  shown 
in  Table  3-2.  A  nine  layer  crust  is  consistent  with  the  crust  assumed 
by  Ellis  and  Basham  (1968)  and  Jensen  (1970)  for  central  Alberta 
spectral  ratios. 

Figure  A-lOa  shows  the  comparison  between  the  theoretical 
spectral  ratios  for  the  seventeen  layer  crust  at  Leduc  and  the  nine 
layer  crust.  They  are  almost  exactly  identical.  The  success  of  this 
experiment  suggests  an  even  further  reduction  of  effective  sedimentary 
layers  to  three  as  used  by  Alpaslan  (1968).  Figure  A-lOb  shows  the 
comparison  between  the  seventeen  layer  model  and  the  three  layer  model 
at  Leduc.  While  there  is  remarkable  similarity  of  the  theoretical 
ratios  at  the  lower  frequencies,  there  are  significant  changes  in  the 
amplitude  of  the  peaks  and  at  the  high  frequencies  the  character  of  the 
ratio  is  quite  different.  This  analysis  resulted  in  the  use  of  the  nine 
layer  model  for  most  efficient  and  reliable  studies  of  the  effect  of 
sediments  on  the  theoretical  ratios. 


Figure  A-ll  shows  the  comparison  between  observed  and  theo- 
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retieal  spectral  ration  based  only  on  the  sedimentary  crust  at  Loduc 
as  modelled  above.  There  is  littie  similarity  between  the  curves. 

This  result  agrees  with  the  findings  of  Ellis  and  Basham  (1968)  who 
explained  this  lack  of  agreement  as  probably  due  to  scattering  at  the 
base  of  the  crust.  Another  possible  explanation,  however,  is  a 
fundamental  error  in  the  crustal  model. 

Analysis  of  the  structure  beneath  the  three  stations  shows 
thickening  of  the  sediments  to  the  southwest.  Fernandez  (1966)  has 
shown  that  all  other  parameters  being  the  same,  a  thickening  of  the 
model  should  cause  the  peaks  of  the  spectral  ratios  to  migrate  toward 
smaller  frequencies  since  those  waves  which  are  able  to  constructively 
and  destructively  interfere  in  the  thicker  crust  must  have  a  longer 
wavelength.  This  effect  is  demonstrated  in  Figure  A-l  for  a  very  simpl 
two  layer  crustal  model.  As  the  thickness  of  this  crust  is  increased 
the  peaks  migrate  to  lower  and  lower  frequencies.  As  predicted  from 


the  thicknesses  of  the  sediments ,  VNG  located  on  the  thickest  sediments 
of  the  three  stations  should  have  peaks  offset  toward  the  lower  frequen 
cies,  OBS  on  the  thinnest  sediments  should  have  peaks  toward  the  higher 
frequencies,  and  OSC  station  should  be  in  between.  Figure  A-12  shows 
the  actual  results  at  the  three  field  stations.  The  various  peaks  do 
not  migrate  as  predicted,  hence  other  factors  than  simple  thickening 
must  be  affecting  the  spectral  ratios  differently  at  the  three  stations 


3.2  Weathering  Layer  Effects  on  the  Spectral  Ratio 

A  possible  explanation  for  the  difference  in  the  observed 
spectral  ratios  at  the  three  stations  and  for  lack  of  agreement  between 


experimental  and  theoretical  results  was  thought  to  be  changes  in  she 
low  velocity  weathering  layer.  Bayrock  and  Hughes  (1962)  have  studie 
the  surficial  geology  of  the  Edmonton  area.  The  bedrock  in  the  Leduc 
area  is  composed  of  upper  Cretaceous  shales,  shaly  sandstones,  and  sa 
stones  which  dip  twenty  feet  per  mile  to  the  southwest.  This  bedrock 
is  overlain  by  eight  to  ten  feet  or  more  of  glacial  drift.  The  three 
stations  are  located  in  an  area  called  the  Gwynne  Outlet  produced  by 
the  erosive  action  of  escaping  waters  from  Pleistocene  Lake  Edmonton. 
A  shaly  sand  deposited  in  the  Gwynne  Outlet  varies  from  one  to  five  f 
in  thickness  and  is  now  on  the  surface  at  the  stations.  The  exact  de 
of  the  till  at  the  stations  was  not  known,  and  a  hammer  seismograph 
survey  was  performed  to  determine  the  depths  and  velocities  of  the 
surficial  material  at  the  OBS  station.  The  upper  ten  feet  was  found 
to  be  very  low  velocity  material,  1250  feet  per  second,  perhaps  the 
shaly  sandy  soil  mentioned  above.  The  interval  from  ten  to  thirty  fe 
in  depth  has  an  average  velocity  of  4300  feet  per  second  and  is  proba 
glacial  till.  6700  feet  per  second  material  was  found  at  a  thirty  fc 
depth  and  is  possibly  bedrock.  The  horizontal  extent  of  this  refract 
survey  was  only  about  four-hundred  feet,  hence,  it  is  possible  that  m 
surficial  layers  than  are  indicated  above  are  present.  Furthermore, 
surficial  layer  of  velocity  lower  than  the  layer  above  it  could  net  h 
been  detected  by  the  refraction  method  used  in  the  survey. 

Spectral  ratios  were  theoretically  computed  for  models  ir.cl 
ing  surficial  layering  and  were  compared  with  the  same  models  without 
the  weathering  layer.  Using  the  hammer  seismograph  results  and  a 
usual  Poisson's  ratio,  almost  no  weathering  layer  effect  was  noted  on 
the  spectral  ratios.  However,  following  Erickson  et  al  (1968),  a  rat 
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of  S  to  P  velocities  of  1/6  was  used  as  an  extreme  case  for  the  uncon¬ 
solidated  surface  layers  to  a  depth  of  .02  kilometers  resulting  in  the 
spectral  ratios  shown  in  Figure  A-13.  Two  observations  can  be  made: 

1)  The  character  of  the  ratio  at  higher  frequencies  is  more 
affected  by  surficial  layering  than  at  lower  frequencies.  This  effect 
was  expected  considering  the  very  short  time  lags  involved  in  multiple 
reflections  within  the  weathering  layer. 

2)  The  effect  of  the  weathering  layer  is  to  lower  the  ampli¬ 
tude  of  the  spectral  ratio  over  most  of  the  frequency  range.  Tnis 
phenomenon  is  mildly  surprising  since  one  would  expect  a  low  velocity 
surface  layer  to  refract  the  incoming  P-coda  into  a  more  vertical 
direction,  hence,  resulting  in  a  higher  vertical  to  horizontal  surface 
motion  ratio.  However,  because  of  the  high  Poisson’s  ratio  possible 

in  unconsolidated  surficial  material ,  it  is  the  converted  S-waves  which 
become  more  important.  When  an  S-wave  is  refracted  toward  the  vertical, 
its  transverse  particle  motion  will  result  in  a  decrease  of  the  ratio  of 
vertical  to  horizontal  surface  motion. 

Tne  weathering  layer  is,  therefore,  seen  as  a  possible  impor¬ 
tant  factor  in  spectral  ratio  studies  in  regard  to  the  amplitude  of  the 
peaks,  but  not  very  significant  as  to  the  relative  position  of  the  peaks. 
Ellis  and  Basham  (1968)  experienced  some  difficulty  in  obtaining  matches 
between  theoretical  and  experimental  ratios  at  several  stations  because 
of  the  high  amplitude  of  their  theoretical  ratios.  It  is  conceivable 
that  the  surficial  layers  at  these  stations  might  have  had  a  consider¬ 
able  effect  on  the  observed  data  and  might  have  altered  spectral  ratio 
amplitudes  considerably. 
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3.3  Poisson's  Ratio  in  Sedimentary  Rocks 

In  previous  work  on  modeling  the  sedimentary  crustal  section 
(Ellis  and  Basham,  1963),  the  rocks  have  been  assumed  to  be  ideal 
Poisson's  solids  with  a  Poisson's  ratio  of  1/4.  There  have  been  few 
actual  field  measurements  to  verify  this  assumption,  however.  If 
Poisson's  ratio  is  not  held  at  1/4,  a  considerable  change  in  the 
theoretical  spectral  ratios  can  be  expected.  Koefford  (1955)  studied 
the  effect  of  Poisson's  ratio  on  the  reflection  coefficients  of  plane 
waves.  His  calculations  for  angles  of  incidence  in  the  range  0°  to 
30°  showed  significant  effects  on  the  reflection  coefficients  which 
increased  with  the  larger  angles  of  incidence  and  decreased  with  smaller 
compressional  velocities.  From  elasticity  theory,  it  is  well  known 
that  the  ratio  of  compressional  to  shear  wave  velocity  in  a  rock  depends 
only  on  Poisson's  ratio.  The  5-wave  velocity  is  seen  to  contribute 
significantly  to  the  spectral  ratio.  Calculations  indicate  that  peaks 
in  the  horizontal  transfer  function  and,  therefore,  troughs  in  the 
theoretical  spectral  ratios  occur  exactly  at  those  frequencies  for 
which  the  S-wave  is  able  to  reinforce  itself  after  one  or  more  reflec¬ 
tions  off  the  free  surface  (Fernandez,  1966).  Hence  Poisson's  ratio 
for  its  effect  both  on  reflection  coefficients  and  on  S-wave  velocity 
is  an  important  factor  to  be  considered  carefully  when  constructing 
crustal  models  for  the  sedimentary  layers. 

Koefford,  Oosterveld,  and  Galons,  (1963),  made  laboratory 
measurements  of  the  P  and  S  velocities  of  twenty-nine  limestones  of 
varying  porosity.  They  found  that  Poisson's  ratio  decreased  with 
increasing  velocity  and  increased  with  water  saturation.  For  water 
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saturated  rocks  of  porosity  less  than  10%,  not  an  untypical  ease  for 
central  Alberta  sediments,  the  average  value  of  Poisson's  ratio  was 
.311.  Z ism an  (1933)  made  laboratory  measurements  on  Poisson's  ratio 
and  concluded  that  "Poisson's  Ratio  increases  greatly  with  mean  stress 
applied",  an  observation  that  could  have  considerable  bearing  on 
assigning  Poisson's  ratios  to  deeply  buried  sedimentary  rocks. 

The  difficulty  in  measuring  Poisson's  ratio  in  situ  is  mostly 
in  efficiently  generating  S -waves  at  the  source.  Oliphant  (1950) 
attempted  to  identify  transverse  body  waves  by  shooting  vertical  pro¬ 
files  with  a  three  component  seismometer  on  the  surface  near  a  borehole 
and  with  shots  at  various  depths  in  the  borehole,  thus,  any  transverse 
motion  on  the  surface  would  have  been  due  to  shear  effects.  The 
experiment  was  not  a  success:  "Evidently  the  fraction  of  detonation 
delivered  to  shear  waves  was  exceedingly  small  in  comparison  with  that 
delivered  by  other  extraneous  motions." 

Richards  (1933)  gives  P  and  S  velocities  from  near  surface 
measurements  in  Southwest  Persia  for  overburdens  consisting  of  sand¬ 
stone,  marls,  anhydrite,  salt,  and  limestone.  The  ratio  of  P  to  S 
velocity  ranged  from  1.8  to  2.3  with  an  average  of  2.0,  which  corre¬ 
sponds  to  a  Poisson's  ratio  of  1/3.  Richards  (1961)  then  proceeded 
to  use  a  Poisson’s  ratio  of  one-third  in  analysing  wide  angle 
reflection  data  in  the  Alberta  foothills. 


Meizner  (1965),  Molotova  (1963),  and  Okado  (1964)  have  also 
made  near  surface  in  situ  measurements  of  Poisson's  ratio.  They 
reported  shear  to  comprcssional  wave  velocity  ratios  as  low  as  1/7 
for  low  velocity  weathering  layers. 

With  the  advent  of  mechanical  vibrators,  significant  progres: 
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has  been  made  in  recent  years  in  determining  S-wave  velocities  at  depth. 
Cherry  and  Waters  (1968)  describe  one  such  experiment,  "In  the  field, 
the  vibrator  generates  horizontal  motion  in  the  ground  while  suitable 
geophones  are  arranged  to  pick  up  this  motion,  SH  waves  are  received 
if  the  source  motion  and  geophones  axes  are  at  right  angles  to  the  spread 
direction.  In  a  simple  parallel  layered  system,  SH  is  reflected  and 
refracted  as  SH  only  and  gives  rise  to  a  simple  reflection  record." 

In  Medford,  Oklahoma  such  shallow  reflections  indicated  a  Poisson's 
ratio  of  .394.  They  found  that  for  most  of  the  reflection  examples, 
the  average  ratio  of  shear  velocity  to  compressional  velocity  in  the 
first  few  thousand  feet  was  near  .5,  that  is,  Poisson’s  ratio  was  equal 
to  1/3. 

The  continuous  signal  method  utilizing  a  mechanical  vibrator 
was  applied  to  borehole  measurements  by  Erickson,  Miller,  and  Waters 
(1968).  They  set  up  the  vibrator  to  generate  S-waves  downward  from  the 
surface  and  lowered  a  geophone  to  various  depths  down  to  ten-thousand 
feet  to  record  the  signal.  Interval  (shear  to  compressional  velocity) 
ratios  were  calculated  from  the  combination  of  downhole  S— wave  data  and 
continuous  velocity  log  data.  The  values  ranged  from  .45  to  .70  with 
an  average  overall  of  .52,  an  average  Poisson's  ratio  slightly  greater 
than  1/3. 

Figure  A-14  shows  the  theoretical  spectral  ratios  for  the 
nine  layer  Leduc  model  for  values  of  Poisson’s  ratio  varying  from  .25 
to  .38.  Superimposed  on  each  of  these  is  the  OBS  field  data.  While 
none  of  the  models  produce  a  satisfactory  fit,  the  models  with  Poisson's 
ratio  about  .33  seem  to  be  best. 

Using  Keofford's  et  al  (1963)  assertion  that  Poisson's  ratio 
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decreases  with  increasing  velocities  as  a  basis,  a  nine  layer  crustal 
model  was  computed  using  a  Poisson’s  ratio  of  .25  in  hard  layers  (P-wave 
velocities  greater  than  4.0  KM/sec.)  and  .33  in  soft  layers  (P-wave 
velocity  less  than  4.0  KM/sec.).  This  model  is  shown  in  Figure  A-15a 
along  with  a  model  all  layers  of  which  have  a  Poisson's  ratio  equal  to 
one-third.  The  difference  between  these  two  models  is  small  indicating 
that  it  is  the  lower  velocity  layers  in  the  sediments  which  are  most 
strongly  influencing  the  spectral  ratios,  at  least,  in  so  far  as  S-wave 
velocity  is  concerned. 

Hence,  it  has  been  shown  that  Poisson's  ratio  can  have  a 
pronounced  effect  on  the  spectral  ratios  and  S-wave  velocities  must  be 
chosen  carefully  when  modelling  the  sedimentary  crust.  However,  it 
must  be  noted  that  no  value  of  Poisson's  ratio  produces  a  match  between 
observed  and  theoretical  ratios  for  models  based  on  only  the  sedimentary 
crustal  layers  at  Leduc. 

3.4  Effect  of  Anisotropy  on  the  Spectral  Ratio 

Another  possible  cause  for  the  lack  of  agreement  betv/een 
theoretical  and  observed  spectral  ratios  might  be  anisotropy  in  the 
crustal  layers.  In  general,  the  array  of  moduli  of  elasticity  relating 
stress  to  strain  in  an  elastic  body  is  given  by  a  symmetrical  six  by 
six  matrix  containing  twenty-one  independent  terms.  The  strain  energy 
function  for  such  a  body  is  a  homogeneous  quadratic  function  of  six 
variables.  These  parameters  may  be  simplified  for  homogeneous  substances 
which  have  some  degree  of  symmetry.  If  the  body  possesses  symmetry  in 
all  directions  so  that  the  strain  energy  function  is  invarient  to  all 
transformations  of  axes, 


the  body  is  said  to  be  isotropic.  Otherwise 
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it  is  anisotropic  or  aleotropic.  A  transversely  isotropic  media 
possesses  an  axis  of  symmetry  in  the  sense  that  all  rays  at  right  angles 
to  this  axis  are  equivalent.  Such  symmetry  is  possessed  geometrically 
by  certain  crystals  of  the  rhombohedral  and  hexagonal  systems.  From 
the  large  scale  seismic  point  of  view,  this  is  the  symmetry  found  in 
flat  sedimentary  layers  and  planar  igneous  bodies.  White  (1955)  and 
Postma  (1955)  have  shown  further  that  a  stratified  section  composed  of 
alternating  layers  of  different  elastic  materials  may  be  considered  as 
a  single  homogeneous  transversely  isotropic  layer  for  wavelengths 
large  compared  to  individual  layer  thicknesses.  Hence,  transverse 
isotropy  can  be  expected  to  arise  in  crustal  layers  from  two  sources, 
either  the  nature  of  the  rock  material  itself  or  the  manner  in  which 
thin  isotropic  rock  layers  are  stratified. 

Stonely  (1949)  was  one  of  the  first  to  study  plane  waves  in 
a  transversely  isotropic  media.  Musgrave  (1954),  Postma  (1955),  and 
Jolly  (1956)  added  the  theory  of  anisotropic  waves  emanating  from  a 
point  source.  These  authors  indicated  that  the  major  differences 
betv/een  wave  propagation  in  transversely  isotropic  media  and  in  purely 
isotropic  media  are  that,  in  the  anisotropic  body,  there  exists  depen¬ 
dence  of  wave  velocity  on  direction  of  travel  and  the  waves  do  not 
separate  into  perfect  P  and  S  types ,  that  is ,  waves  for  which  the 
divergence  or  curl  of  the  displacement  vanish  respectably.  A  derivation 
of  velocities  in  a  transversely  isotropic  medium  is  given  in  the  appendix. 

The  fact  that  rocks  are  not  isotropic  has  been  observed  for 
many  years.  Koefford  et  al  (1963)  determined  longitudinal  and  transverse 
velocities  of  twenty-nine  limestones  of  varying  porosity  and  the  elastic 
anisotropy  factor,  that  is,  the  ratio  of  transverse  velocity  to  longi- 
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tudinal  velocity,  was  found  to  vary  over  a  wide  range  from  0.5 
a  third  of  which  were  over  1.05.  Layat  et  al  (1961)  noted  that  in  order 
to  satisfactorily  identify  refractions  from  marker  beds  in  the  Sahara, 
it  was  necessary  to  take  anisotropy  into  account,  for  without  anisotropy 
depth  calculations  put  marker  beds  in  the  middle  of  known  homogeneous 
beds  without  velocity  differentiations.  Cherry  et  al  (1968)  from  stucie 
of  shear  wave  propagations,  determined  an  anisotropic  factor  of  1.08 
for  the  upper  sedimentary  section  in  Medford,  Oklahoma.  Dunoyer  de 
Segonzac  and  Leherrere  (1959)  determined  anisotropic  factors  in  the 
northern  Sahara  ranging  from  1.00  in  sandstones  to  1.20  in  andhydrite 
by  a  method  utilizing  a  velocity  log  together  with  the  first  arrival 
times  recorded  at  several  geophone  levels  from  shotpoints  at  various 
distances  from  the  well.  Vanderstoep  (1966)  utilized  a  similar  method 
in  New  Mexico  and  Texas.  The  red  beds  in  a  2400  foot  well  in  Lea  County 
New  Mexico  showed  anisotropic  factor  of  1.07,  the  salt  beds  1.00.  On 
the  Texas  Gulf  coast,  a  12,000  foot  well  yields  anisotropic  factors  of 
1.03  or  less.  McCollum  (1932),  Weatherby  (1934),  Ricker  (1953),  and 
White  and  Sengbus’n  (1953)  working  in  Canada,  Oklahoma,  Colorado,  and 
Texas  respectively,  observed  near  surface  anisotropy  in  shales  and 
limestones  ranging  from  1.1  to  1.4.  J.  Cholet  and  H.  Richard  (1954) 
determined  an  anisotropic  factor  of  1.09  in  the  northern  Sahara  down 
to  1250  meters.  A.H.  Kleyn  (1956)  employed  the  results  from  a  refrac¬ 
tion  profile,  a  well  velocity  survey,  and  a  radial  well  survey  in 
Betun,  South  Sumatra.  He  concluded  that  the  anisotropic  factor  was 
1.00  at  the  surface,  1.15  at  1500  meters,  and  1.04  deeper.  He  further 
noted  that,  at  great  depths,  the  horizontal  velocity  may  locally 
decrease  or  remain  constant  even  though  the  vertical  velocity  increases 


with  depth.  Anderson  (1961)  noted  that  anisotropy  ray  be  respar, 
for  certain  discrepancies  in  surface  wave  data  including  disagre 
between  Love  wave  and  Rayleigh  wave  studies  of  both  velocity  arc 
tide  motion. 


sible 

par- 


In  order  to  calculate  the  theoretical  spectral  ratio  for  bod 
waves  in  a  system  of  transversely  isotropic  layers,  the  matrix  method 
of  Haskell  (1953)  can  be  used  provided  the  terms  involving  elastic 
constants  are  properly  modified  beforehand.  Anderson  (1961)  and  Ander 
and  Harkrider  (1962)  were  the  first  to  apply  the  Haskell  method  to 
transversely  isotropic  media  for  the  purpose  of  doing  surface  v:ave 
dispersion  studies.  A  brief  outline  of  the  application  of  the  Haskell 
matrices  to  body  waves  in  a  transversely  isotropic  layered  system  as 
conceived  by  E.R.  Kanasewich  of  The  University  of  Alberta  is  given  in 
the  appendix.  This  derivation  leads  to  a  matrix  formula  for  calculati 
the  theoretical  ratio  for  P-waves  which  is  particularly  suitable  for 
the  digital  computer. 

Due  to  the  low  angles  of  incidence  which  P-waves  from  tele- 
seismic  events  impinge  on  the  crust,  it  is  reasonable  to  expect  that 
the  effect  of  anisotropy  on  the  spectral  ratios  will  be  small  for  such 
events  since  the  direction  of  wave  propagation  is  essentially  vertical 
in  these  cases,  and  the  horizontal  velocities  are  only  slightly  trough 
into  play.  Figure  A-lSa  shows  the  comparison  of  the  theoretical  spect 
ratios  for  an  isotropic  model  and  an  anisotropic  model  for  the  sedimen 
tary  layers  at  Leduc  for  an  anisotropic  factor  of  1.2  and  an  angle  of 
incidence  of  20°.  While  a  fair  amount  of  change  in  amplitude  of  the 
peaks  is  noted,  there  is  almost  no  frequency  shifting  of  the  peaks. 

the  same  models  at  an  angle  of  incidence  of  40^ . 


Figure  A-16b  shows 
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Once  again  a  significant  change  in  the  amplitude  of  the  peaks  is  evident, 
but  not  so  much  as  to  change  the  overall  character  of  the  ratio.  How¬ 
ever,  at  an  angle  of  incidence  of  60°  as  shown  in  Figure  A- 16c,  both 
the  amplitude  and  position  of  the  peaks  are  severely  affected  by 
anisotropy. 

At  larger  angles  of  incidence,  however,  anisotropy  is  impor¬ 
tant  and  must  be  taken  into  consideration.  Several  researchers 
(Alpaslin,  Ellis  and  Basham)  have  noted  difficulties  in  obtaining 
reasonable  matches  between  theoretical  and  experimental  ratios  for  near 
events.  Formerly  this  has  been  explained  by  anomalous  SV  conversions 
and  simultaneous  arrivals  of  energy  from  different  ray  paths.  To  these 
theories  must  be  added  the  strong  effect  of  anisotropy  on  weaves  arriving 
at  such  large  angles  of  incidence. 

Figure  A-15b  shows  the  comparison  of  the  isotropic  and  aniso¬ 
tropic  models  at  Leduc  for  an  angle  of  incidence  of  26°  which  is  thought 
to  be  the  actual  case  for  earthquake  epicenters  in  the  Aleutians,  The 
difference  is,  as  expected,  mostly  only  a  matter  of  amplitude  although 
at  a  frequency  of  just  under  2.0  Hz  a  small  peak  does  appear  on  the 
anisotropic  ratio  which  was  not  present  on  the  isotropic  ratio.  It  is 
evident  that  at  this  angle  of  incidence  it  is  not  likely  that  anisotropy 
could  be  the  cause  of  any  large  discrepancies  in  character  between 
theoretical  and  observed  spectral  ratios. 

3,5  Effect  of  Deep  Crustal  Layering 

It  is  probable  that  the  sedimentary  structure  alone  can  account 
for  neither  the  general  character  of  the  observed  spectral  ratios  nor 


the  change  in  the  character  of  the  ratios  fro:r.  station  to  station. 
Before  disregarding  the  spectral  ratio  method,  however,  one  should 
consider  the  possibility  of  deeper  layers  in  the  crust  affecting  the 
spectral  ratios.  Ellis  and  Basham  (1968)  have  shown  that  the  additie 
of  the  "gross"  crustal  layers,  that  is,  those  layers  between  the  sedi 
ments  and  the  Mohorovicic  discontinuity,  does  not  have  a  significant 
effect  on  the  general  characteristics  of  the  spectral  ratios.  The 
deep  crustal  layers  added  were  a  generalization  of  those  given  by 
Kanasewich  and  Cumming  (1966)  as  the  results  of  regional  refraction 
studies  in  southern  Alberta.  Figure  A-17a  shows  the  theoretical  spec 
tral  ratios  for  a  model  consisting  only  of  the  sediments  at  Leduc  and 
for  a  model  of  the  sediments  plus  the  deeper  crustal  layers  used  by 
Ellis  and  Basham.  The  two  ratios  are  quite  similar  in  general  char¬ 
acteristics.  The  effect  of  the  deep  layers  is  the  addition  of  small 
oscillations  on  the  general  curve  shape.  These  small  oscillations 
are  not  apparent  on  a  similar  comparison  given  by  Ellis  and  Basham 
(1968)  probably  because  of  the  large  amount  of  smoothing  performed  in 
their  analysis.  The  presence  of  the  small  oscillations  is  in  agree¬ 
ment  with  the  work  of  Phinney  (1964)  who  predicted  that  any  layer 
would  produce  oscillations  on  the  spectral  ratios  with  frequency 
spacing  equal  to  one  over  the  associated  time  lag. 

Hence  the  deeper  crustal  layers  because  of  their  long  time 
lags  (the  two  way  time  to  the  Mohorovicic  discontinuity,  for  example, 
is  about  twelve  seconds)  should  produce  oscillations  with  very  short 
frequency  spacings.  This  apparent  failure  of  the  deeper  crustal 
layers  to  affect  the  spectral  ratio  coupled  with  the  lack  of  corre¬ 
lation  success  with  spectral  ratios  based  on  only  the  sediments,  leav 
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only  one  possibility  for  successful  application  of  the  method,  that  is, 
intermediate  crustal  layers  in  the  upper  few  kilometers  of  the  pre- 
Cambrian  rocks. 

3.6  Models  Involving  a  Shallow  Thin  Soft  Layer  in  the  Pre-Cambrian 

The  pre-Cambrian  in  Alberta  has  been  studied  by  Burwash  (1957), 
Garland  and  Burwash  (1959),  and  Kanasewich  and  Clowes  (1969).  Burwash 
(1957)  presented  data  from  petrological  studies  of  pre-Cambrian  core 
samples  showing  that  the  rock  types  in  the  basement  vary  considerably 
in  different  localities.  A  well  in  the  Leduc  area  actually  contained 
two  different  rock  types  in  the  upper  few  feet  of  the  basement  complex. 
Garland  and  Burwash  (1959)  combined  petrological  and  geophysical  data 
to  produce  a  rough  lithological  map  of  the  basement  in  central  Alberta. 
Kanasewich  and  Clowes  (1959),  utilizing  deep  seismic  sounding  and 
gravity  data  in  southern  Alberta  mapped  a  rift  in  the  pre-Cambrian 
extending  down  to  the  Mohorovicic  discontinuity.  The  pre-Cambrian  in 
Alberta  is,  hence,  seen  to  be  a  complex  geological  system  rather  than 
a  homogeneous  crystalline  substance  as  often  supposed.  It  is  hot 
unreasonable,  therefore,  to  suggest  the  existance  of  layering  in  the 
uppermost  basement  rocks  which  could  have  a  significant  effect  on  the 
spectral  ratios. 

An  initial  attempt  to  arrive  at  a  suitable  theoretical  spectral 
ratio  to  match  the  observed  data  at  Leduc  by  means  of  an  additional 
layer  In  the  basement  was  based  on  three  factors:  1)  Visual  inspection 
of  the  observed  ratios  indicated  that  a  better  match  could  be  obtained 
if  the  theoretical  peaks  were  shifted  toward  lower  frequencies. 

Fernandez  (1965)  noted  that  it  is  possible  to  perform  such  a  shift  by 


increasing  the  depth  of  the  model.  One  method  of  increasing  the 
is  the  addition  of  another  layer  on  the  bottom  of  the  sediments. 

2)  The  gravity  map  of  central  Alberta  shows  the  Leduc  area  to  be 
located  in  the  center  of  a  broad  gravity  low.  A  shallow  "soft”  1 
in  the  basement  could  explain  such  a  low.  3)  Commercial  seismic 
reflection  operations  in  the  Edmonton  area  often  fail  to  produce 
clear  reflection  from  the  basement  surface.  A  possible  explanati 
this  phenomenon  could  be  a  low  reflection  coefficient  at  the  top 
''soft"  basement  layer. 

Several  hundred  theoretical  spectral  ratios  were  general 
using  models  with  an  additional  layer  on  the  bottom  of  the  sedimer 
the  best  of  which  is  shown  in  Figure  A-17b  along  with  the  observe: 
ratio  at  OSC  station.  The  fit  is  apparently  reasonable.  The  mode 
involves  a  soft  layer  of  about  2000  feet  thickness  beneath  the  sec 
Unfortunately,  no  suitable  models  were  found  to  fit  the  data  for  t 
other  stations  as  well,  hence,  the  fit  for  the  OSC  data  is  probabl 
nothing  more  than  a  good  illustration  of  the  ambiguity  of  the  inve 
process. 

3.7  A  Three  Second  Marker  in  the  Pre-Cambrian 

Evidence  for  a  possible  marker  bed  in  the  pre-Cambrian  a 
two-way  time  of  3.0  to  3.5  seconds  comes  from  two  independent  sour 
Referring  back  to  Figure  A- 2 ,  the  cepstrum  analysis  of  observed  sp 
tral  ratios  at  OSC  stations,  one  can  note  an  interesting  phenomena 
The  graphs  show  the  effect  which  difference  in  time  lags  have  on  t 
spectral  ratios.  The  time  lag  which  would  probably  have  the  most 


on  the  observed  ratio  is  the  two-way  reflection  time  for  various  markers 
in  the  layered  system.  The  three  graphs,  each  for  a  different  event 
recorded  at  the  same  stations,  show  remarkable  similarities  for  short 
time  lags,  that  is,  short  two-way  times  undoubtably  associated  with 
sedimentary  layers.  It  is  very  interesting  that  the  similarity  between 
the  three  events  continues  out  to  about  3.5  seconds,  while  the  total 
two-way  time  in  the  sediments  around  Leduc  is  less  than  1.6  seconds. 

One  explanation  for  the  significant  time  lag  at  3.5  seconds  could  be 
multiples  within  the  sediments,  but  a  more  interesting  possibility  is 
the  existance  of  a  true  marker  in  the  pre-Cambrian  at  a  two-way  time 
of  around  3.5  seconds. 

Recent  deep  seismic  reflection  records  from  the  area  several 
miles  north  of  the  Leduc  Observatory  obtained  by  G.  Gumming  of  The 
University  of  Alberta  have  indeed  shown  a  distinct  reflection  at  a 
two-way  time  of  3.4  seconds. 

The  generation  of  theoretical  spectral  ratios  for  models  with 
a  marker  in  the  3.0  to  3.5  second  range  can  lead  to  various  possibilities 
for  the  thickness  and  P— velocity  of  the  basement  layer  since  it  is  only 
the  ratio  of  the  terms,  the  transit  time,  which  is  known.  Hence,  the 
spectral  ratio  method  could  be  of  considerable  use  in  determining  the 
actual  thickness  of  the  basement  layer,  if  the  various  models  do  not 
prove  to  be  completely  non-unique. 

At  the  Leduc  Observatory,  no  satisfactory  fit  was  obtained 
v/ith  any  model  involving  a  3.4  second  marker,  however,  a  reasonable  fit 
was  obtained  with  a  3.1  second  marker  as  shown  in  Figure  A- 18a.  The 
essential  feature  of  this  model  is  the  low  P-wave  velocity  used  in  the 
basement  layer  (4.6  KM/sec.).  Tnis  low  velocity  keeps  the  reflection 


coefficient  at  the  bottom  of  the  sediments  small.  Figure  A-19b  shows 
the  theoretical  ratios  for  a  two-way  time  of  3.1  seconds  for  models 
with  higher  P-velocities  in  the  basement.  The  models  become  more  and 
more  like  the  theoretical  ratio  for  a  model  based  only  on  the  sediments 
as  the  reflection  coefficient  at  the  basement  surface  is  increased. 

Further  encouragement  for  the  3.0  second  marker  models  comes 
from  a  reasonable  fit  for  the  VNG  station  as  shown  in  Figure  A-19a. 

While  the  fit  is  by  no  means  excellent,  it  is  by  far  the  closest  of 
any  of  the  hundreds  of  models  in  this  study  attempted  for  the  VNG 
station.  This  particular  model  is  for  a  pre-Cambrian  marker  at  3.0 
seconds  or  only  0.1  second  different  from  the  OBS  station  which  is 
located  some  twrelve  miles  away. 

To  see  if  the  lower  crustal  layers  might  improve  the  fit  for 
the  3.1  seconds  model  at  OBS,  the  presumed  Conrad  and  Mohorovicic 
discontinuities  we re  added  to  the  model ,  and  the  resultant  spectral 
ratios  are  shown  in  Figures  A-18b  and  A-18c.  The  closely  spaced 
modulations  due  to  these  layers  do  not  alter  the  theoretical  spectral 
ratios  to  an  important  degree. 

No  fit  was  found  for  the  OSC  station,  however,  which  is 
located  between  the  other  two  stations,  six  miles  distant  from  each. 

This  lack  of  a  match  may  indicate  a  great  deal  of  complicated  structure 
in  the  pre-Cambrian  basement  in  the  Leduc  area,  or  it  could  also  point 
out  the  inherent  non-uniqueness  and  consequent  lack  of  certainty  in  the 
matches  at  the  OBS  and  VNG  stations.  A  detailed  gravity  survey  is 
currently  being  undertaken  in  hopes  of  gaining  more  information  about 
the  sub-basement  structures,  if  any,  in  the  Leduc  area. 

The  models  which  most  closely  fit  the  OBS  and  VNG  data  produce 
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theoretical  spectral  ratios  which  are  too  high  in  amplitude.  Ellis  and 
Basham  (1968)  have  encountered  a  similar  problem,  tentatively  explaining 
the  low  observed  ratios  as  due  to  scattering  and  anomalous  P  to  S 
conversions  near  the  base  of  the  sediments.  Another  possible  explana¬ 
tion  might  be  as  discussed  in  section  3.2,  effects  of  non— consolidated 
near  surface  sediments  with  an  extremely  high  Poisson's  ratio.  At  any 
rate,  the  problem  of  the  amplitude  difference  between  observed  and 
theoretical  spectral  ratios  in  the  Leduc  area  remains  unsolved  at  this 
time. 


3.8  Future  Spectral  Ratio  Studies 

The  spectral  ratio  method  in  the  high  frequency  range  has  been 
found  to  be  potentially  very  sensitive  to  two  parts  of  the  crust  which 
are  often  overlooked  in  other  types  of  seismic  investigations,  namely, 
weakly  consolidated  surficial  layers  and  relatively  shallow  structures 
in  the  crystalline  basement.  The  method  could  prove  to  be  a  very 
valuable  mapping  tool  for  these  types  of  studies  in  the  future. 

For  investigations  of  the  lower  crust,  however,  the  high 
frequency  spectral  ratio  inversion  is  unsuitable.  Future  spectral  ratio 
studies  of  the  deep  crust  should  probably  involve  longer  period  data 
from  only  carefully  selected  events  in  order  that  the  effects  of  shallow 
crustal  layering  and  of  pP  type  phases  might  be  minimized. 

One  of  the  more  interesting  lessons  to  arise  from  these  studies 
into  the  spectral  ratio  method  was  the  surprising  difficulty  encountered 
in  obtaining  even  a  few  plausible  matches  between  theoretical  and 
observed  data  considering  the  limited  frequency  range,  the  considerable 
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number  of  free  parameters,  and  the  inherent  non-uniqueness  of  the 
problem.  'This  observation  indicates  that  the  most  efficient  and 
instructive  method  of  handling  such  inversions  in  the  future  might  be 
based  on  a  computerized  curve  matching  technique  since  the  apparent 
scarcity  of  good  fits  would  make  any  good  match,  based  on  a  model 
consistent  with  the  known  geology,  worth  further  investigation. 
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FIGURES 
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Figure  A-l.  Theoretical  spectral  ratios  for  simple  two  layer 

models  at  various  thicknesses  (D)  but  with  identical 
layer  velocities.  The  thicker  the  model,  the  closer 
the  spectral  peaks  approach  one  another. 
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Figure  A-2.  Time  lag-spectral  ratio  curve  for  the  OSC  events. 

The  0SC3  and  0SC4  spectral  ratios  show  peaks  at 
about  12  seconds  which  corresponds  to  their  pP 
phase  arrival  time.  0SC2  with  a  pP  time  of  6 
seconds  does  not  show  the  distinctive  peaks  at 


12  seconds. 
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Figure  A-3.  Comparison  of  the  spectral  ratios  for  the  Oct.  31 
event;  one  with  a  ten  second  P-coda  window,  the 
other  with  a  forty  second  window.  The  forty  second 


window  provides  considerably  more  resolution. 
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Figures  A-4,  A-5,  A-6.  Observed  spectral  ratios  at  the  OBS ,  VNG, 

and  OSC  stations.  The  extent  to  which  the 
observed  ratios  are  alike  at  each  station 
is  a  measure  of  the  reliability  of  the 


spectral  ratio  method 
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Figure  A- 7. 


Comparison  of  1969  data  at  OBS  with  data  collected 


by  other  researchers  in  1965  and  1968  also  at  OBS. 
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Figure  A-8. 


Map  showing  the  area  of  The  University  of  Alberta’ 


Leduc  Observatory. 
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Fi giare  A-9.  Theoretical  spectral  ratios  for  two  nine  layer 
sedimentary  models,  the  only  difference  between 
them  being  the  density  distribution.  Trie  density 
is  seen  to  affect  only  the  height  of  the  peaks, 


not  their  frequency. 
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Figure  A-lOa.  The  comparison  between  the  spectral  ratios  for  a 

seventeen  layer  model  at  Leduc  and  a  nine  layer 
model.  They  are  almost  identical. 

Figure  A-lOb.  The  comparison  between  the  seventeen  layer  model 

and  a  three  layer  model.  There  are  significant 


differences  between  the  two. 


Figure  A— 11.  Comparison  between  observed  and  theoretical  spectral 
ratios  using  a  model  based  only  on  the  sedimentary 
crustal  layers. 


seikauJ 


to 


SEDIMENTARY 

MODEL 


0 


1 


Figure  A— 12.  Comparison  of  the  observed  spectral  ratios  at  the 
three  field  stations.  The  various  peaks  do  not 
migrate  as  predicted  from  the  known  geology,  hence 
other  factors  than  simple  sedimentary  thickening 
must  be  affecting  the  spectral  ratios  differently 


at  the  three  stations. 


Figure  A-13. 


Theoretical  spectral  ratios  for  models  with  and 


without  the  weathering  layer. 
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Figure  A-14.  Theoretical  spectral  ratios  for  the  nine  layer 
Leduc  model  for  values  of  Poisson's  ratio 
varying  from.  0.25  to  0.38.  Superimposed  on 


each  of  these  is  the  OBS  data. 
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Figure  A— 15a*  Theoretical  spectral  ratio  for  a  sedimentary 

crustal  model  with  a  Poisson's  ratio  of  0.25 
in  hard  layers  and  0.33  in  soft  layers  compared 
with  the  spectral  ratio  for  the  same  model  with 
Poisson's  ratio  equal  to  0.33  in  all  the  layers. 

Figure  A-15b.  Comparison  of  spectral  ratios  for  the  isotropic 

and  anisotropic  sedimentary  crustal  models  at 


Leduc . 
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Figure  A- 16a,  A-16b,  A-16c.  Comparison  of  spectral  ratios  for 

isotropic  and  anisotropic  sedimentary 
crustal  models  at  various  angles  of 
incidence. 
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Figure  A-17a.  Theoretical  spectral  ratios  for  a  model  consisting 

of  only  the  sediments  at  Leduc  and  for  a  model 
consisting  of  the  sediments  plus  the  deeper  crustal 
layers. 

Figure  A-17b.  Theoretical  and  observed  spectral  ratios  at  station 

OSC  for  a  model  involving  a  shallow  2000  foot  layer 


beneath  the  sediments. 
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Figure  A— 18.  Comparison  of  the  observed  ratio  at  OBS  with  models 
involving  a  3.1  second  marker,  a  3.1  second  marker 
plus  the  Conrad  discontinuity,  and  a  3.1  second 
marker  plus  both  the  Conrad  and  Mohorovicic 
discontinuities . 


Figure  A- 19a.  Observed  spectral  ratio  at  VNG  station  compared 

with  the  theoretical  spectral  ratio  for  a  model 
involving  a  3.0  second  marker  in  the  basement. 


Figure  A- 19b.  Theoretical  spectral  ratios  for  models  with  a 

3.1  second  marker  but  with  varying  basement 


velocity. 
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APPENDIX  B 
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wi 


IN  A  TRANSVERSELY  ISOTROPIC  MEDIA 


B.l  Introduction 


A  mathematical  discussion  of  seismic  wave  propagation  in  a 
transversely  isotropic  media  and  the  derivation  of  the  matrix  formula 
for  the  P-coda  spectral  ratio  in  a  horizontally  layered  anisotropic 
system  are  given  in  this  appendix.  Listed  below  are  all  the  math¬ 
ematical  symbols  used  in  order  of  appearance  as  well  as  a  qualitative 
description  of  each. 


x,  y,  z,  or  x± 

x *  ,  y* ,  z ' 

0 

c.  . 
ij 

e.  . 
ij 

W 

P.  . 

u.  or  u,  v,  w 
1 

u.  or  u,  v,  w 
l 

t 


r . 

l 

k. 

i 


Cartesian  coordinates 

Transformed  coordinates 

Angle  of  rotation 

Matrix  of  elastic  constants 

Strain  tensor 

Strain  energy  function 
Stress  tensor 

Displacement  in  the  i  direction 
Particle  velocity  in  the  i  direction 
Time 

Direction  cosines  for  type  i  waves 
Cotangent  of  angle  of  incidence 
Wave  number  for  type  i  waves 
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v_  . 

1 

Velocity  for  type  i  waves 

A,  y 

Lame ' s  Constants 

c 

Phase  velocity  along  surface 

p 

Density 

u,,  V1A 

d  d 

Maximum  amplitudes  of  downgoing  wave  in 

the  horizontal  and  vertical  directions 

03 

Angular  frequency 

a  ,  a 
v  H 

P-wave  velocity  in  vertical  and 

horizontal  directions 

3  ,  3U 

v  H 

S-wave  velocity  in  vertical  and 

horizontal  directions 

UdCt’  Ud$’  ^uO' ’  Uu3  Maximum  amplitude  of  the  downgoing  and 


upgoing  "P"  and  "S"  waves  respectively 

k 

Horizontal  wave  number 

P,  Q,  R,  S 

Exponential  factors  for  sinusoidal  plane 

waves 

*V  iY3 

Ratios  of  vertical  and  horizontal  maximum 

amplitudes  for  "P"  and  "S"  waves  respec¬ 
tively;  functions  only  of  the  elastic 

parameters,  density,  frequency,  and 

angle  of  incidence 

U  U  U  U4  Sums  and  differences  of  maximum  amplitudes 


X  ,  XQ,  Y  ,  YQ 
a’  3  a’  3 

of  upgoing  and  downgoing  "P"  and  "S" 

waves 

Stress  amplitude  factors;  functions  only 

of  the  elastic  parameters,  density, 

frequency,  and  angle  of  incidence 
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m,  m-1, 

m-2,  etc.  Subscripts  indicating  layer  under 

analysis  as  counted  from  the  top  of  the 

layered  system 

D 

m 

4x4  matrix  relating  the  particle 

velocity  and  stress  vector  to  the 

amplitude  vector  in  layer  m 

V 

m 

Particle  velocity  and  stress  vector  for 

layer  m 

E 

m 

Matrix  D  with  z  =  0 
m 

E  -1 
m 

Inverse  of  matrix  E 

m 

A  ,  B 
nr  m 

Combinations  of  stress  and  displacement 

amplitude  factors:  A  is  always 

m 

imaginary,  is  always  real 

a 

m 

-1 

Product  of  matrices  D  and  E  for 

m  rn 

layer  m 

n 

Number  of  layers  in  the  system 

J 

Matrix  product  for  all  layers  combined 

U 

m 

Amplitude  vector  for  layer  m 

D 

Determinant  of  the  coefficients  from  the 

equations  relating  J  to  the  boundary 

conditions  of  the  media 

R 

P-coda  spectral  ratio 

U  ,  W 

P  P 

Horizontal  and  vertical  transfer  functions 

B.2  Equations  of  Motion 


For  a  transversely  isotropic  medium  with  a  vertical  (z)  axis, 
the  following  directional  and  rotational  substitutions  should  not 
affect  the  nature  of  seismic  propagation;  that  is,  they  should  not 
alter  the  strain  energy  function,  W, 


x  *  -  x 

y  ’  =  y 

z'  =  -z  B.l 


due  to  symmetry  about  any  horizontal  plane,  and 

x  =  x’  cos  0  -  y’  sin  0 
y  =  x’  sin  0  +  y’  cos  0 
z  =  z1 

B.2 


due  to  symmetry  about  the  vertical  axis  where  the  primed  symbols  are 
the  transformed  coordinates  and  0  is  the  angle  of  rotation.  Love 
(1920)  performed  these  substitutions  in  the  strain  energy  function  and 
determined  that  twenty-one  independent  elastic  constants  are,  in  fact, 
reduced  to  five  in  a  transversely  isotropic  body.  That  is, 
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c  .  . 
13 


c 


c  c  0  0  0 


1 1 

1  2 

1  3 

1  2 

C 

1  1 

c 

1  3 

0 

0 

0 

1  3 

c 

1  3 

c 

3  3 

0 

0 

0 

0 

0 

c 

44 

0 

0 

0 

0 

0 

c 

4  4 

0 

0 

0 

0 

0 

c 
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where  c 


c  -  c 
11  12 


6  6 


The  strain  energy  function  for  a  transversely  isotropic 
medium  reduces  to  the  following: 
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B.4 


where  e. .  is  the  strain  tensor, 
ij 

The  stresses  are  then: 
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The  equations  of  motion  are  given  by : 
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3  P  .  . 
13 


3  x . 


B.3  Derivation  of  Velocities  in  a  Transversely  Isotropic  Medium 

In  order  to  simplify  the  problem  to  a  two  dimensional  ca 
motion  will  only  be  allowed  in  the  x-z  plane.  SH  waves  are,  there 
excluded,  the  displacement  in  the  y  direction  is  zero,  all  derivat 
with  respect  to  y  are  zero,  and  the  equations  of  motion  become: 


3 2  u 
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xy 


+ 


0  X 
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32w  9P 

p  -  = - 5* 

3  t2  3x 

To  introduce  the  anisotropic  condition,  equations  B.5  are 
substituted  in  the  equations  of  motion  as  follows: 
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zz 
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32u 
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3  3u 

-  (  c  -  +  c 
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3  z 
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3z 
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3  f 


3x 


{  (  c  ( 
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du 

3z 
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3w 

3x 


5 }  +  ~sr  ( C13 
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3x 


+  c 


3w 


3  3 
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3z 
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Suitable  trial  solutions  for  downgoing  plane  waves  with  harmonic 
time  dependence  will  be  of  the  form 


u  -  Ud  •  exp  (  iwt  -  i  k.  (l±x  +  ruz  )  ) 

w  =  -W  •  exp  (  iwt  -  ik.  (  l.x+n.z  )  )  B.9 

d  ill 

where  1.  and  n.  and  k.  are  the  direction  cosines  and  the  wavenumber 
ill 

for  v/aves  of  type  i. 

Substituting  this  trial  solution  B.9  in  the  equations  of 
motion,  B.8,  gives  a  new  equation  of  motion  which  may  be  written  in 


matrix  form  as 


0 


(l.2c 


+  n .  c 

1  11  1  44 


pc.,2} 


1  .  II  .  ( c  +  c  ) 

11  44  13 


-l.n.  (c  +  c  ) 

11  44  13 


-  {l.2c  +  n. 2c  -  pc. 2} 

144  133  1 


U 


where  cw  is  the  velocity  for  waves  of  type  i. 

The  velocities  of  wave  propagation  can  be  determined  by 

2 

solving  for  the  two  roots  in  c^  by  setting  the  determinant  of  tne 
coefficients  in  equation  B.10  equal  to  zero. 

pc . 2}  -  1 . 2n. 2 (c  +  c  ) 0 

1  11  13  44 

B.  11 

For  vertical  rays  1.  =  0,  n.  =  1,  or 

l  l 


{1/ 
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-t  n 


1 1 


4  4 


pc.2} 

1 


(l.: 

1 


+  n. 


44 


3  3 


a 

V 


c 

-  3  3—  =  Vertical  P  velocity 

P 


c 

-  -4-4—  =  Vertical  S  velocity 

P 


For  horizontal  rays,  1  =  1,  n^  =  0,  or 


c 

1 1 


P 


Horizontal  P  velocity 


c 

o  4  4 

P  -  — -  =  Horizontal  SV  velocity 

V  P 


and  in  a  three  dimensional  case : 


c  ~  c 

=  - — - - =  Horizontal  SH  velocity  B.  12 

p 


For  a  wave  in  any  arbitrary  direction 


c . 
1 


2 


1 


2p 


c  1 . 2  +  c  n  .  2  +  c 
-111  331  41* 


C 

4  4 


2 


c  )  n.2  }2  +  4  (  c  +  c  )2  l.V2!  1/2 

44  1  13  44  1  3.  J 
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Both  solutions  of  this  equation  are  real  and  positive.  If  the  three 
dimensional  case  were  being  considered,  another  wave  would  exist  as 
follows,  shown  by  Stonely  (1949). 


)  +  n . 2c  } 
1  44 


B.  14 


Hence  there  are  three  different  velocities  to  be  found  for  a  wave 
propagating  in  any  direction. 

To  identify  these  waves,  it  is  useful  to  consider  the  special 
case  of  the  purely  isotropic  medium  in  which 
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c  =  c  =  A  +  2u 

11  3  3 


c  =  c 
13  12 


C 

4  4 


P 


where  A  and  y  are  Lame's  Constants. 

The  solutions  of  equation  B.13  and  B. 14  then  yield: 


B.  15 


A  +  2y 

B.  16 

B.  17 

B.  18 

Thus  in  a  transversely  isotropic  body,  SV  and  SH  waves 
generally  travel  at  different  velocities. 


or  a  P  wave 


y 

P 


or  a  SV  wave 


or  a  SH  wave 


B.4  Theoretical  Spectral  Ratios  for  P-Waves  in  a  System  of  Horizontal 
Transversely  Isotropic  Layers 


A  solution  to  the  elastic  wave  equation  valid  for  a  transvesely 
isotropic  layer  containing  both  upward  and  downgoing  waves  is  as  following 


=  {U,  e 

da 


-iP 


TT  lPl 

U  e  } 
ua 


R  .  / 
e  +  1 


U  e  ^  +  U 

d3 


u3 
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7c. 


=  -lY 


f  II 
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da 


-ir  lP  I  R  .  r  TT  -iQ  TT  iQ 

-  -  U  e  }  e  -  iyQ  {  U,Qe  -  U  e 

ua  p  dp  ua 


where:  U,  ,  'J,  .  U  ,  U  represent  the  maximum  amplitudes  of  the 

da’  dg  ua’  u^  F  * 

downgoing  and  the  upgoing  "P"  and  "S"  waves  respectively. 

P,  Q,  R,  S  provide  the  plane  wave  characteristics  as  follows: 

P  =  kr  z 

a 

Q  =  kr^z 

R  =  i  (  wt 
S  =  i  (  wt 

iYa » iy ^  are  the  ratios  of  vertical  and  horizontal  amplitudes  for  P 
and  S  waves  respectively  as  can  be  derived  from  equations  B.10. 

For  waves  of  type  i : 
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a  a 
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di 


U,. 

di 
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c  +  c  r .  -  pc  r.(c  +  c  ) 

11  i(i)  1  1  13  44 


-  r.  (  c  +  c  )  -  {c  +c  r.2-pc2} 

1  13  44  44  331 


The  following  substitutions  are  necessary  to  convey  the  solution  into 
a  more  useful  form: 
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Since  the  position  on  the  surface  is  arbitrary  for  plane  waves,  x  can 
be  set  to  zero.  The  quantity  measured  is  not  displacement  but  particle 
velocity,  so  the  derivative  of  equation  B.19  yields 


— =  ik  {  U  cos  P  +  iU  sin  P  +  U  cos  Q  +  iU  sin  Q  } 

C  2  14  3 

— —  =  ik  {-  U  y  sin  P  +  ill  Y  cos  P  -  U  yD  sin  Q  +  iU  yQ  cos  Q}  _ 

c  2  a  i  a  4p  3P  B.  2  x 

Combining  equations  B.19,  B.20,  and  B.5  leads  to  stress  relationships 
as  follows : 


P  =  ik  {  iX  U  cos  P  -  X  U  sin  P  +  iXnU  cos  Q  -  XQU  sin  Q} 

zz  a  2  a  i  P  4  p  3 

P  =  ik  {  iY  U  sin  P  +  Y  U  cos  P  +  iYQU  sin  Q  4-  YGU  cos  Q} 

ZX  a  2  P4  b3 
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X  ,  X  are  purely  imaginary  while  y  Yn  are  real, 
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Equations  B.21  and  B„22  can  be  expressed  in  matrix  form  for  any  layer  m. 
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where  is  the  amplitude  vector,  and  V  is  the  particle  velocity  and 


stress  vector. 
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=  0,  and  to  preserve  continuity 
:  interface,  the  following  relation 
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3.23  and  B.24,  one  has 
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where 
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Then  B.21  becomes  for  the  (m-1)  layer: 
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Repeatedly  substituting  in  equation  B.20,  a  recurrent  relationship  can 
be  established  starting  with  the  bottom  layer  n  until  the  surface  is 
reached : 
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Using  substitutions  B.20,  the  boundary  condition  that  stress  must 
vanish  at  the  free  surface,  and  assuming  an  incident  P  wave  at  the 
base  of  layer  n-1 
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Solving  these  equations  for  the  normalised  particle  velocities, 
one  arrives  at  the  crustal  transfer  functions. 
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vertical  transfer  function 


where  D  is  the  determinant  of  coefficients  from  equations  B.31 


D  —  —  (J  +  J  )  (J  +  J  )  +  — —  (J  +  J  )  (J  +  J  ) 

k2  12  22  31  41  k2  11  21  32  42 


The  ratio  of  the  moduli  of  the  above  complex  expressions  gives  the 
theoretical  P-coda  spectral  ratio  (R)  for  any  assumed  crustal  model. 

It  can  be  shown  that  elements  j  and  J  are  purely  Imaginary  whereas 
d  and  J  are  real. 
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The  spectral  ratio  can  thus  be  seen  to  be  dependent  only  on  the  matrix 
J,  a  function  of  the  layer  parameters,  frequency,  and  angle  of  incidence 


of  the  wave. 


APPENDIX  C 


P-CODA  PLOTS  OF  THE  POUR  FALL  1969  ALEUTIAN  EVENTS 
AS  RECORDED  AT  STATIONS  OBS ,  VNG ,  AND  OSC 
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